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摘要 科学 研究 中 的 计算 机 模拟 称 为 科学 模拟 (scientific simulation) ， 文 章 从 其 狭义 和 广义 层面 出 发 ， 把 科 
学 模拟 分 为 数值 计算 (numerical computation), #447 4& (simulation intelligence) 和 科学 大 脑 (science 
brain) 3 个 阶段 ， 并 描述 了 各 阶段 的 特征 。 目 前 ， 科 学 模拟 正在 进入 模拟 智能 阶段 ， 即 在 科学 大 数据 和 人 工 
智能 的 驱动 下 ， 科 学 模拟 正 从 传统 的 数值 计算 逐步 转向 与 人 工 物 能 相 融 合 的 模拟 方式 。 文 章 针对 支撑 模拟 物 


能 阶段 的 计算 系统 展开 讨论 ， 阅 述 了 其 设计 指导 思想 、 


基本 方法 和 关键 技术 等 问题 。 
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1 科学 研究 中 的 计算 机 模拟 


人 类 进入 大 科学 时 代 后 ,，“ 模 拟 ” 作 为 “理论 ” 
“实验 ”以 外 重要 的 补充 技术 手段 ， 成 为 科学 研究 的 
第 3 个 支柱 。 从 表现 形式 看 ， 科 学 研究 可 以 被 视 为 建 
立 模型 (modeling) 的 过 程 。 而 模拟 (simulation) 则 
是 所 建立 的 科学 模型 在 计算 机 上 的 运行 过 程 。 最 早期 
的 计算 机 模拟 (computer simulation) 可 以 追溯 到 第 二 
次 世界 大 战 之 后 ， 是 专门 针对 核 物理 学 和 气象 学 研究 
的 一 种 开创 性 的 科学 工具 。 后 来 ， 在 越 来 越 多 的 学 科 
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中 ， 计 算 机 模拟 变 得 愈 发 重要 ， 并 不 断 衍 生出 计算 和 
其 他 领域 交叉 的 学 科 ， 如 : 计算 物理 、 计 算 化 学 和 计 
算 生 物 学 等 学 科 。WeaverD 在 1948 年 撰文 指出 : 人 类 
解决 有 序 复杂 问题 并 实现 科学 新 飞跃 ， 将 主要 依赖 于 
计算 机 技术 的 发 展 和 不 同学 科 背 景 科 学 家 的 技术 碰 
撞 。 一 方面 ， 计 算 机 技术 的 发 展 使 人 类 能 够 解决 复杂 
且 难 以 处 理 的 问题 。 另 一 方面 ， 计 算 机 技术 能 够 有 效 
激发 有 序 复 杂 性 问题 的 新 解决 方法 。 这 种 新 解决 方法 
本 身 也 正 是 计算 科学 (computational science) 的 范畴 
之 一 ， 让 科学 家 有 机 会 集中 资源 ， 将 不 同 领域 的 见解 
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聚焦 在 共同 问题 上 。 这 种 见解 聚焦 的 结果 ， 促 使 不 同 
学 科 背 景 的 科学 家 们 ， 形 成 比 单一 学 科 背 景 科学 家 们 
更 强大 的 “混合 团队 ”; 这 样 的 “混合 团队 ”将 有 能 
力 解决 某 些 复杂 性 问题 ， 并 且 得 出 有 用 的 结论 。 总 
之 ， 科 学 和 建 模 紧 密 相 关 ， 模 拟 执行 代表 理论 的 模 
型 ， 人 们 把 科学 研究 中 的 计算 机 模拟 称 为 科学 模拟 
(scientific simulation) 。 

目前 ， 还 没有 任何 针对 “计算 机 模拟 ”的 单一 定 
义 能 够 恰当 地 描述 科学 模拟 的 概念 。 美 国 国防 部 将 模 
拟定 义 为 一 种 方法 ， 即 :“ 一 种 随时 间 实 现 模型 的 方 
iE". 进而 ， 将 计算 机 模拟 定义 为 一 种 过 程 ， 即 :“ 在 
计算 机 上 执行 代码 、 控 制 和 显示 接口 硬件 ， 并 与 现实 
世界 设备 进行 接口 交互 的 过 程 ”。WinsbergP 把 计算 机 
模拟 的 定义 又 分 为 狭义 和 广义 范围 。 

在 狭义 定义 中 ， 计 算 机 模拟 就 是 “在 计算 机 上 运 
行程 序 的 过 程 ”。 计 算 机 模拟 使 用 步 进 方法 来 探索 数 
学 模型 的 近似 行为 。 模 拟 程序 在 计算 机 上 的 一 次 运行 
过 程 ， 代 表 了 对 目标 系统 的 一 次 模拟 。 人 们 愿意 用 计 
算 机 模拟 方法 求解 问题 ， 主 要 有 以 下 2 个 原因 : © 原 
始 模型 本 身 包含 离散 方程 ; @) 原始 模型 的 演化 更 适合 
用 “规则 ”， 而 不 是 “方程 ”来 进行 描述 。 值 得 注意 
的 是 ， 这 种 狭义 角度 指 代 计 算 机 模拟 时 ， 需 要 限定 到 
特定 处 理 器 硬件 上 算法 的 实现 、 用 特定 编程 语言 编写 
应 用 ,以 及 核 函数 程序 、 使 用 特定 编译 絮 等 限制 条 
件 。 在 不 同 应 用 问题 的 场景 下 ， 由 于 这 些 限制 条 件 的 
变化 ， 通 常会 获得 不 同 的 性 能 结果 。 

在 广义 定义 中 ， 可 以 把 计算 机 模拟 看 作 研 究 系 统 
的 一 种 综合 方法 ， 是 更 加 完整 的 计算 过 程 。 该 过 程 包 
括 模 型 选择 、 通 过 模型 实现 、 算 法 输出 计算 、 结 
数据 可 视 化 及 研究 。 整 个 模拟 的 过 程 也 可 以 与 科学 
研究 过 程 进行 对 应 ， 如 LynchB 所 描述 : D 提出 一 个 
经 验 上 可 回答 的 问题 ; @ 从 旨 在 回答 该 问题 的 理论 
中 推导 出 一 个 可 证 伪 的 假设 ; Co 收集 (或 发 现 ) 和 
分 析 经 验 数据 以 检验 该 假设 ; 0D 拒绝 或 未 能 拒绝 该 
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假设 ; © 将 分 析 结 果 与 得 出 该 问题 的 理论 联系 起 来 。 
在 过 去 ， 这 种 广义 的 计算 机 模拟 通常 出 现在 认识 论 或 
者 方法 论 的 理论 场景 中 。 

Winsbere" 进一步 将 计算 机 模拟 划分 为 基于 方程 
的 模拟 (equation-based simulation) 和 基于 主体 的 模 
拟 (agent-based simulation)。 基 于 方程 的 模拟 常用 于 
物理 等 理论 学 科 中 。 这 些 学 科 中 一 般 存在 主导 性 的 理 
论 ， 这 些 理论 可 以 用 来 指导 构建 基于 微分 方程 的 数学 
模型 。 例 如 ， 基 于 方程 的 模拟 可 以 是 针对 粒子 的 模 
拟 ， 这 种 模拟 通常 包含 数量 巨大 的 多 个 独立 粒子 和 一 
组 描述 粒子 之 间 相 互 作用 的 微分 方程 。 此 外 ， 基 于 方 
程 的 模拟 也 可 以 是 基于 场 的 模拟 ， 通 常 包含 一 组 描述 
连续 介质 或 场 的 时 间 演 化 方程 。 基 于 主体 的 模拟 往往 
遵循 某 种 演化 规则 ， 是 模拟 社会 和 行为 科学 的 最 常见 
Z AX. fü. Schelling 的 隔离 政策 模型 中 。 尽 管 基于 
主体 的 模拟 在 一 定 程度 上 可 以 表示 多 个 主体 的 行为 ， 
但 与 基于 方程 的 粒子 模拟 不 同 ， 这 里 没有 控制 粒子 运 
动 的 全 局 微分 方程 。 

从 计算 机 模拟 的 定义 和 分 类 中 ， 可 以 看 出 人 们 对 
科学 模拟 不 同 层 次 的 期 望 。 从 狭义 的 计算 机 模拟 角度 
看 ， 它 已 经 成 为 理论 分 析 和 实验 观察 等 传统 认 知 方式 
的 补充 手段 。 科 学 或 工程 领域 无 一 例外 是 由 计算 机 模 
拟 推动 的 ， 在 某 些 特定 应 用 领域 和 场景 下 ， 甚 至 是 由 
计算 机 模拟 改变 的 。 如 果 没 有 计算 机 模拟 ， 许 多 关键 
技术 就 无 法 被 理解 、 开 发 和 利用 。 广 义 的 计算 机 模拟 
蕴含 着 一 个 哲学 问题 ,计算 机 是 否 可 以 自主 进行 科学 
研究 ? 科学 研究 的 目标 是 认识 世界 ， 这 意味 着 计算 机 
程序 必须 创造 新 的 知识 。 随 着 人 工 智能 技术 研究 及 应 
用 的 新 一 轮 爆 发 ， 人 们 对 计算 机 自动 地 以 “智能 ” 方 
式 进 行 科学 研究 充满 了 期 待 。 值 得 一 提 的 是 ，KitanoDl 
在 2021 年 提出 的 “ 诺 贝 尔 -图 灵 挑 战 ” 的 新 观点 一 一 
“到 2050 年 ， 开 发 能 够 自主 执行 研究 任务 的 智能 科学 
家 ， 做 出 诺 贝 尔 奖 级 别 的 重大 科学 发 现 ”。 该 观点 涉 
及 狭义 和 广义 的 计算 机 模拟 相关 技术 ,但 没有 围绕 广 


义 定义 的 “哲学 问题 ”深入 探讨 ， 只 是 把 其 作为 科学 
模拟 的 一 个 宏伟 目标 看 待 。 


2 科学 模拟 的 发 展 阶 段 


从 最 直观 的 视角 来 看 ， 科 学 模拟 的 载体 是 计算 机 
程序 。 从 数学 形式 上 讲 ， 计 算 机 程序 是 由 可 计算 函数 
组 成 的 ， 其 中 每 个 函数 将 有 限 输入 数据 的 离散 集 映 射 
到 有 限 输出 数据 的 离散 集 上 。 从 计算 机 技术 上 讲 ， 计 
算 机 程序 等 于 算法 加 上 数据 结构 。 因 此 ， 科 学 模拟 的 
实现 需要 以 科学 问题 及 其 解决 方式 被 形式 化 抽象 为 条 
件 。 这 里 ， 本 文 借用 Simon 的 观点 : 科学 家 即 问题 
“求解 器 "。 在 此 观点 中 ， 科 学 家 给 自己 设 定 了 重大 科 
学 问题 ， 确 定 问题 和 解决 问题 的 策略 和 技术 是 科学 发 
现 的 本 质 。 基 于 上 述 “ 求 解 器 ”的 话语 体系 ， 本 文 类 
比 求解 方程 的 形式 ， 将 科学 模拟 的 发 展 划分 为 3 个 阶 
段 ， 即 数值 计算 、 模 拟 智能 和 科学 大 脑 〈 图 1)。 

2.1 数值 计算 

计算 机 模拟 最 初 用 于 执行 科学 研究 过 程 所 需要 数 
值 计 算 。 因 此 ， 在 许多 口头 或 正式 的 学 术 文章 中 ， 人 
们 往往 把 数值 计算 、 数 值 模拟 和 计算 机 模拟 不 加 区 分 
地 统一 进行 使 用 和 表述 。 第 一 台电 子 数字 计算 机 
ENIAC 是 用 来 计算 弹道 方程 ， 进 而 用 于 风 洞 设计 和 天 
气 预报 。 在 传统 的 数值 计算 阶段 ， 计 算 机 模拟 程序 执 
行 的 过 程 就 是 基于 数理 方程 的 求解 过 程 ， 即 简单 的 
SREE y= F(x)。 事 实 上 ， 在 过 去 的 数 十 年 中 ， 科 
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学 研究 领域 的 数值 计算 推动 了 超级 计算 机 运算 速度 的 
不 断 提升 。 人 们 通过 人 硬件、 软件 和 算法 的 创新 来 持续 
提高 y = F(x) 的 求解 速度 和 精度 。 正 因 如 此 ， 通 常 是 
以 求解 稠密 或 者 稀 跑 线性 代数 方程 的 测试 程序 ， 作 为 
评估 超级 计算 机 速度 的 基准 测试 程序 ， 如 Linpack 中 和 
HPCGM 基 准 测试 程序 。 

然而 ， 这 种 将 部 分 复杂 科学 问题 转换 为 相对 简单 
的 计算 问题 的 解 题 模 式 ， 仅 仅 是 一 种 粗 粒 度 的 建 模 方 
案 ， 在 一 些 应 用 场景 下 会 遇 到 计算 瓶 贷 。 在 解决 真实 
场景 中 复杂 物理 模型 时 ， 常 常 面临 基本 物理 原理 计算 
量 过 大 的 问题 ， 并 由 此 导致 空 有 原理 而 无 法 有 效 解决 
科学 问题 。 例 如 ， 第 一 性 原理 分 子 动力 学 的 关键 为 求 
解 量子 力学 Kohn-Sham 方 程 ， 其 核心 算法 求解 过 程 是 
多 次 求解 大 规模 特征 值 问 题 。 特 征 值 问 题 的 计算 复杂 
ENN (N 为 矩阵 的 维度 )。 在 实际 物理 问题 的 求解 
中 ， 最 常用 的 平面 波 基 组 通常 是 原子 个 数 的 100 一 
10 000 倍 。 这 意味 着 对 于 上 千 原 子 的 体系 规模 ， 和 矩阵 
维度 N 达 到 10"， 其 相应 的 浮 点 数 运算 总 量 也 将 达到 
10 FLOPS， 即 达到 EFLOPS 级 别 的 计算 量 ?。 需 要 注 
意 的 是 ， 在 单 步 分 子 动力 学 中 需要 多 次 求解 特征 值 问 
题 ， 这 也 就 使 得 单 步 分 子 动力 学 的 模拟 时 间 通 常 为 数 
分 钟 乃 至 1 小 时 。 由 于 单 步 分 子 动力 学 的 模拟 物理 时 
间 只 能 达到 1 飞 秒 ,假设 要 完成 纳 秒 物理 时 间 的 分 子 
动力 学 模拟 过 程 ， 就 需要 10 个 分 子 动力 学 步 。 相 应 
的 计算 量 至 少 要 达到 10” FLOPS。 如 此 庞大 的 运算 量 


I 用 物理 原理 和 海量 数据 构建 智能 模型 端 到 端的 智能 科学 活动 构建 新 理论 


数值 计算 


(numerical computation ) 


y= Ff(x),a(x)) 


模拟 智能 


(simulation intelligence) 


科学 大 脑 


(science brain ) 


图 1 科学 模拟 发 展 的 3 个 阶段 
Figure 1 Three phases of scientific simulation 
F/If 表 示 机 理 建 模 ，A/a 表示 基于 人 工 智能 建 模 
F/f. mechanistic modeling, A/a: AI-based modeling 


(D FLOPS ,每 秒 浮 点 运算 次 数 ;EFLOPS ,每 秒 进 行 百 亿 亿 次 浮 点 运算 。 
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即使 使 用 世界 上 最 大 规模 的 超级 计算 机 也 难以 在 短 时 
间 内 完成 。 为 了 解决 仅 使 用 第 一 性 原理 计算 带 来 的 超 
大 计算 量 , 研究 人 员 发 展 了 多 尺度 方法 ， 其 中 最 典型 
的 是 获得 了 2013 年 诺 贝尔 化 学 奖 的 量子 力学 /分 子 力 
学 (QM/MM) 方法 。 该 方法 的 思想 是 针对 核心 物理 
化 学 反应 部 分 (如 : 酶 及 其 结合 底 物 的 活跃 位 点 原 
T), 采用 高 精度 的 第 一 性 原理 计算 方法 ， 对 于 周围 
的 物理 化 学 反应 区 域 (溶液 、 蛋 白质 和 其 他 区 域 ) 采 
用 低 精 度 且 计算 复杂 度 更 低 的 经 典 力学 方法 。 这 种 高 
精度 、 低 精度 相 结 合 的 计算 方法 ， 可 以 有 效 地 降低 计 
量 。 但 面 对 实际 问题 时 ， 该 方法 依然 存在 着 巨大 的 
挑战 。 例 如 ， 细 胞 半径 约 0.2 微米 的 单个 生殖 支原体 
包含 3x10? 个 原子 和 77 000 个 蛋白 质 分 子 。 由 于 核心 
计算 时 间 仍 来 自 QM 部 分 ， 模 拟 2 小 时 的 过 程 预计 需 
要 耗费 10" 年 。 如 果 将 类 似 问 题 推 广 到 人 脑 的 模拟 中 ， 
相应 的 系统 原子 数 将 达到 105 个 ， 保 守 估 计 需 要 101 
个 活跃 位 点 的 QM 计算 。 由 此 可 以 推断， 模拟 1 小 时 
的 QM 部 分 需要 长 达 10” 年 的 时 间 ， 而 MM 部 分 的 模 
拟 也 需要 长 达 10” 年 的 时 间 。 这 种 超 长 计算 时 间 的 情 
况 也 被 称 为 “维度 灾难 ”。 
2.2 模拟 智能 

如 前 所 述 ， 传 统 数值 计算 在 探索 微观 世界 的 多 尺 
度 模拟 上 遭遇 “维度 灾难 ”"， 即 随 着 维 数 的 增加 ， 计 
算 代价 呈 指 数 级 增加 。 著 名 物理 学 家 Dirac" 指出 : 大 
部 分 物理 学 和 全 部 化 学 的 数学 理论 所 需 的 基本 物理 定 
律 是 完全 已 知 的 。 但是， 这 些 问 题 的 困难 在 于 这 些 定 
律 的 精确 求解 方法 包含 过 于 复杂 的 方程 求解 过 程 ， 以 
至 无 法 在 有 限时 间 内 完成 求解 。 巴 斯 德 实验 室 研究 团 
BUH 2022 年 撰写 的 调研 报告 中 提出 了 “模拟 智能 ” 
(simulation intelligence) 的 概念 ， 即 科学 计算 和 人 工 
智能 的 融合 ， 旨 在 研究 以 计算 机 模拟 系统 演化 来 更 好 
地 理解 和 发 现 原 位 现象 。 这 里 的 关键 挑战 是 对 “ 原 
位 ”的 探索 。 因 此 ， 相 比 第 一 阶段 的 科学 模拟 ， 模 拟 
智能 必须 拓展 传统 数值 计算 的 能 力 边 界 ， 发 展 革新 技 
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R (如 融合 人 工 智 能 ) 以 应 对 物理 建 模 和 数据 建 模 中 
的 维度 灾难 。 这 种 对 “ 原 位 ”的 探索 问题 ， 也 正 是 鄂 
维 南 院士 团队 所 倡导 的 科学 智能 (AI4S) "要 解决 的 
关键 问题 。 图 1 中 ， 模 拟 智 能 是 在 传统 的 数值 计算 过 
程 中 庶 入 人 工 智 能 计算 过 程 4(X)， 表 示 成 求解 器 y= 
(f(x),a(x))， 其 显著 特点 是 利用 物理 原理 和 海量 
数据 来 构建 人 工 智能 模型 及 人 工 智能 计算 过 程 。 

因此 ， 模 拟 智能 在 传统 数值 计算 中 炭 人 人工 智 能 
模型 (当前 主要 是 深度 学 习 模 型 ) ， 不 同 于 其 他 人 工 
智能 应 用 领域 的 深度 学 习 模 型 “ 黑 盒子 ”。 模 拟 智能 
要 求 这 些 模型 的 基本 出 发 点 和 基本 结构 是 可 解释 
的 中 。 目 前 ， 这 一 方向 已 存在 大 量 人 研究 ，Zhang 等 中 
在 2023 年 对 模拟 智能 领域 最 新 进展 进行 了 系统 性 的 梳 
理 。 从 理解 亚 原子 ( 波 函 数 和 电子 密度 )、 原 子 (分 
子 、 蛋 白质、 材料 和 相互 作用 ) 到 宏观 (流体 、 气 候 
和 地 下 ) 尺度 物理 世界 ， 把 研究 对 象 分 为 量子 
E F 和 连续 介质 
(continuum) ZKR, MWETA., WEZA, 
小 分 子 、 蛋 白质 、 材 料 科学 、 分 子 间 相互 作用 和 连续 
力学 等 7 个 科学 领域 。 此 外 ， 还 详细 讨论 了 其 中 关键 
的 共同 挑战 ， 即 : 如 何 通 过 深度 学 习 方 法 捕 提 物理 第 
一 性 原理 ， 特 别 是 自然 系统 中 的 对 称 性 。 利 用 物理 原 
理 的 智能 模型 几乎 已 经 渗透 了 传统 科学 计算 的 所 有 和 领 
域 。 模 拟 智能 大 幅 提升 了 对 微观 多 尺度 系统 的 模拟 能 
力 ， 为 在 线 实验 反馈 迭代 提供 了 更 加 全 面 的 支撑 条 
件 。 例 如 ， 计 算 模拟 系统 和 机 器 人 科学 家 之 间 的 快速 
实时 迭代 ， 有 助 于 提升 科研 效率 。 因 此 ， 模 拟 智能 在 
一 定 程 度 上 ， 还 将 包括 “理论 一 实验 ”迭代 的 控制 过 
程 ， 同 时 也 会 涉及 部 分 广义 的 科学 模拟 。 
2.3 科学 大 脑 

传统 的 科学 方法 从 根本 上 塑造 了 人 类 探索 自然 和 
科学 发 现 的 分 步 “ 指 南 ”。 面 对 全 新 的 研究 问题 ， 科 
学 家 们 已 经 被 训练 成 从 假设 和 替代 方案 的 角度 出 发 ， 
指定 如 何 开 展 控制 测试 5 的 定 势 思维 。 虽 然 这 种 研究 
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过 程 在 过 去 几 个 世纪 内 都 是 有 效 的 ， 但 是 非常 缓慢 
的 。 从 某 种 意义 上 来 说 ， 这 种 研究 过 程 是 主观 的 ， 是 
由 科学 家 的 聪明 才智 和 偏见 驱动 的 。 这 种 偏见 ， 有 时 
会 阻碍 必要 的 范式 转变 09。 人 工 智 能 技术 的 发 展 激发 
了 人 们 对 科学 和 智能 融合 产生 最 优 的 且 具 有 创新 性 的 
解决 方案 的 期 望 。 

这 一 阶段 是 对 广义 计算 机 模拟 的 高 级 实现 ， 依 赖 
于 通用 人 工 智能 技术 的 发 展 。 如 图 1 所 示 ， 这 种 高 级 
实现 表示 为 y= 4(f (x))， 意味 着 机 器 (智能 体 ) 能 
够 执行 “创造 、 发 现 ” 新 理论 的 过 程 。 本 文 不 讨论 其 
具体 的 技术 特征 ， 只 描述 宏观 愿景 。Kitano 外 在 提出 类 
似 的 智能 科学 家 的 理念 时 ， 以 2 个 获得 诺 贝尔 奖 的 工 
作 一 一 体 细胞 重 编 程 技术 和 导电 聚合 物 的 发 现 为 例 。 
在 一 定 程度 上 ， 这 2 项 科学 发 现 的 过 程 等 同 于 1 项 相 
同 的 研究 发 现 过 程 ， 即 :“ 有 一 个 明确 的 目标 ， 具 有 
明显 的 科学 意义 ,然后 进行 了 搜索 和 优化 ”"。 尤 其 是 
导电 聚合 物 的 发 现 更 具有 戏剧 性 和 偶然 性 ， 这 一 发 现 
本 质 上 来 源 于 一 场 实验 事故 : Shikawaka 实验 室 的 一 
名 实习 研究 员 错误 地 使 用 异常 高 浓度 的 化 学 物质 来 形 
成 薄膜 ，Shirakawa 注 意 到 这 个 偶然 的 发 现 并 优化 了 薄 
膜 形 成 的 条 件 。 最 终 Shirakawa, MacDiarmiad 和 
Heeger 确 定 了 导电 聚合 物 形 成 的 条 件 。 由 此 可 知 ， 在 
通用 人 工 智能 逐步 成 熟 的 条 件 下 ， 如 果 能 制定 科学 发 
现 的 过 程 ， 并 开发 一 个 可 扩展 的 系统 来 执行 这 个 过 
程 ， 可 能 会 得 到 另 一 种 血 新 的 科学 发 现形 式 ， 即 “ 科 
学 大 脑 ”。 这 种 科学 发 现形 式 将 实现 计划 、 推 测 、 实 
验 、 确 认 和 分 析 等 的 智能 全 自动 化 流程 。 

以 上 所 提 到 的 科学 模拟 发 展 经 历 的 3 个 阶段 ， 能 
够 明显 区 分 计算 机 模拟 在 可 计算 和 智能 化 能 力 方面 逐 
步 提升 的 过 程 。GD 数值 计算 阶段 ， 对 复杂 科学 问题 中 
相对 简单 的 计算 问题 进行 了 粗 粒 度 建 模 ， 属 于 单纯 的 
狭义 计算 机 模拟 定义 范畴 。 它 不 仅 促进 众多 领域 宏观 
尺度 科学 发 现 ， 同 时 也 开启 了 对 微观 世界 的 初步 探 
Ro QD 模拟 智能 阶段 ， 将 针对 微观 世界 的 多 尺度 探索 
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推 上 一 个 新 的 人 台阶。 除了 在 狭义 计算 机 模拟 定义 范畴 
内 对 计算 能 力 进 行 了 数量 级 地 提升 之 外 ， 该 阶段 还 涉 
及 对 实验 中 某 些 关键 环节 的 计算 加 速 ， 在 一 定 程 度 上 
为 科学 模拟 下 一 阶段 的 实现 莫 定 了 基础 。(3) 科学 大 脑 
阶段 ， 将 是 对 广义 计算 机 模拟 定义 的 实现 。 在 此 阶段 
中 ， 计 算 机 模拟 将 具备 创造 知识 的 能 


3 设计 模拟 智能 计算 系统 的 关键 问题 


按照 本 文 对 科学 模拟 发 展 阶段 的 粗 粒度 划分 ， 与 
其 相应 的 计算 系统 也 在 同步 进化 。 超 级 计算 机 在 数值 
计算 阶段 发 挥 了 不 可 替代 的 作用 ;发 展 到 新 的 模拟 智 
能 阶段 ， 底 层 计算 系统 的 设计 也 是 基石 。 那 么 ， 模 拟 
智能 计算 系统 的 发 展 方向 应 该 遵循 什么 样 的 指导 
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纵 观 计算 和 科学 人 研究 发 展 历史 ， 可 归纳 出 计算 系 
统 发 展 的 基本 周期 性 规律 : 在 新 的 计算 模式 和 需求 产 
生 阶 段 的 初期 ， 计 算 系 统 的 设计 侧重 追求 极致 的 专用 
性 。 而 在 经 过 一 段 时 间 的 技术 演变 和 应 用 拓展 之 后 ， 
计算 系统 的 设计 开始 侧重 于 对 通用 性 的 追求 。 在 人 类 
科技 文明 早期 发 展 的 漫长 过 程 中 ,计算 系统 曾经 是 各 
式 各 样 的 专用 机 械 设备 ， 辅 助 进行 一 些 简 单 的 运算 
(图 2)。 近 代 以 来 ， 电子 技 术 的 突破 催生 了 电子 计算 
机 的 出 现 ， 并 且 随 着 其 计算 能 力 的 不 断 提升 ， 数 学 、 
物理 等 学 科 的 发 展 也 不 断 向 前 ， 尤 其 是 超级 计算 机 上 
的 大 规模 数值 模拟 成 果 ， 引 领 了 大 量 前 沿 科 学 研究 和 
重大 工程 应 用 。 由 此 可 见 ， 日 益 发 展 的 通用 高 性 能 计 
算 机 在 不 断 地 加 速 宏观 尺度 科学 的 各 类 大 规模 应 用 ， 
并 取得 重大 成 果 。 接 下 来 ， 微 观 世 界 的 多 尺度 探索 将 
是 未 来 Z 级 〈102) 超级 计算 机 应 用 的 核心 场景 。 而 现 
有 通用 高 性 能 计算 机 的 技术 路 线 则 将 遇 到 功 耗 和 效率 
等 瓶颈 ， 难 以 为 继 。 结 合 模拟 智能 阶段 所 呈现 的 新 特 
征 ， 本 文 认 为 面向 模拟 智能 的 计算 系统 ， 将 以 追求 极 
致 的 Z 级 计算 专用 智能 系统 为 设计 目标 ， 未 来 性 能 最 
高 的 计算 系统 将 专门 针对 模拟 智能 应 用 程序 ， 在 硬件 
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专题 : 大 力 推进 科研 范式 变革 


本 身 及 软件 底层 的 算法 和 抽象 中 进行 定制 。 

直观 上 讲 ， 面 向 模拟 智能 的 计算 系统 离 不 开 智能 
组 件 (软件 和 硬件 )， 那 么 基于 现 有 的 智能 组 件 来 构 
建 智能 计算 系统 就 能 真正 满足 模拟 智能 的 需求 吗 ? 答 
案 是 否定 的 。 李 国 杰 院 士 曾 经 指出 :“ 有 人 曾 戏 让 目 
前 信息 领域 的 形势 为 :“ 软 件 在 吞噬 世界 ， 人 工 智 能 
TETEWEAAE, VRBECECITETRUR A LE, GPU (图 形 
处 理 器 ) ETRE ” BEBE UEBER 
GPU 或 类 似 的 硬件 加 速 器 ， 似 乎 成 了 对 付 大 数据 的 主 
要 出 路 。 但 是 如 果 不 清楚 该 在 什么 地 方 加 速 ， 只 盲目 
依靠 硬件 的 蛮 力 是 不 明智 的 。 因 此 ， 设 计 智能 系统 的 
关键 在 于 深刻 理解 要 求解 的 问题 。 计 算 机 架构 师 的 角 
色 是 选择 好 的 知识 表示 、 识 别 开 销 密集 型 任务 、 学 习 
元 知识 、 确 定 基本 操作 后 ， 再 用 软 硬 件 优化 技术 去 支 
持 这 些 任务 。” 

面向 模拟 智能 的 计算 系统 设计 是 一 个 新 产生 的 研 
究 主题 ， 相 对 其 他 计算 系统 设计 而 言 ， 更 加 具有 显著 
的 独特 性 。 因 此 ， 需 要 一 个 整体 统一 的 视角 ， 来 推进 
人 工 智 能 和 模拟 科学 的 交叉 。1989 年 ，Wah 和 Li05 总 
结 了 关于 智能 计算 机 系统 设计 的 3 个 层次 ， 该 观点 至 
今 依然 值得 借鉴 。 但 遗憾 的 是 ， 目 前 还 没有 任何 关于 
这 方面 的 更 加 深入 的 讨论 和 实际 性 研究 。 具 体 而 言 ， 
& 表示 层 


智能 计算 机 系统 的 设计 要 考虑 3 个 层次 
(representation level)、 控 制 层 (control level) 和 处 理 层 
(processor level) 。 表 示 层 处 理 用 于 解决 给 定 人 工 智 能 
问题 的 知识 和 方法 ， 以 及 如 何 表示 该 问题 ， 控 制 层 关 


宏观 


注 算 法 中 依赖 关系 和 并 行 性 ， 以 及 问题 的 程序 表示 ; 
处 理 层 解决 执行 算法 和 程序 表示 所 需 的 硬件 和 体系 结 
构 组 件 。 下 面 将 以 这 3 个 层次 为 基础 ， 讨 论 面向 模拟 
智能 的 计算 系统 设计 的 关键 问题 。 

3.1 表示 层 

表示 层 是 设计 过 程 中 的 一 个 重要 元 素 ， 包 括 领域 
知识 表示 和 共性 特征 (元 知识 ) 表示 ， 其 决定 了 给 定 
问题 是 否 能 够 在 合理 的 时 间 内 得 到 解决 。 定 义 表示 层 
的 本 质 是 对 适应 广泛 应 用 的 行为 和 方法 进行 高 级 抽 
象 ， 将 它们 与 特定 的 实现 解 厢 。 下 面 给 出 领域 知识 表 
示 和 共性 特征 表示 的 案例 。 

从 现 阶段 面向 科学 的 人 工 智能 研究 看 ， 对 称 性 的 
研究 将 成 为 表征 学 习 的 一 个 重要 突破 口 ， 其 原因 在 于 
物理 上 的 守恒 定律 是 由 对 称 性 导致 的 ( 诺 特 定理 )， 
而 守恒 定律 党 被 用 来 研究 粒子 的 基本 属性 和 粒子 之 间 
的 相互 作用 。 物 理 上 的 对 称 性 是 指 在 某 种 变换 后 或 某 
种 操作 下 的 不 变性 ， 无 法 做 出 可 辨别 的 测量 (不 可 区 
分 性 )。 基 于 多 层 感知 机 (MLP) 、 卷 积 神经 网 络 
(CNN)、 图 神经 网 络 (GNN) 的 小 分 子 表 征 模型 在 有 
效 加 入 对 称 性 之 后 ,已 经 广泛 应 用 于 和 蛋白质、 分 子 、 
晶体 等 物质 的 结构 预测 中 。 

2004 4E, Colella! 向 美国 国防 高 级 研究 计划 局 
(DARPA) 提出 了 科学 计算 的 “七 个 小 矮人 ”(Seven 
Dwarfs) MERRE ARARE Ag 
格 计算 、 非 结构 网 格 计算 、 谱 方法 、 粒 子 方法 、 蒙 特 
卡 治 模拟 。 其 中 ， 每 一 种 科学 计算 问题 ， 都 代表 了 一 
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R 
经 典 
物理 模型 


电子 计算 机 
繁荣 时 代 


通用 电子 计算 机 


高 性 能 计算 
繁荣 时 代 


通用 高 性 能 计算 机 


图 2 科学 模拟 计算 系统 发 展 的 周期 性 规律 
Figure2 Periodic trends of computing system for scientific simulation 
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种 可 以 捕获 计算 和 数据 移动 模式 的 计算 方法 。 受 此 启 
发 ， 巴 斯 德 实验 室 的 Lavin 等 以 类 似 方式 定义 了 模 
拟 智能 中 的 9 种 基 元 (nine motifs of simulation 
intelligence) 一 一 多 物理 现象 多 尺度 建 模 、 代 理 建 模 
仿真 、 基 于 模拟 的 推理 、 因 果 关 系 建 模 推理 、 基 于 主 
体 的 建 模 、 概 率 编程 、 微 分 编程 、 开 放 式 优化 、 机 器 
编程 。 这 9 种 基 元 代表 了 互 为 补充 的 不 同 计算 方法 类 
型 ， 为 协同 模拟 和 人 工 智 能 技术 促进 科学 发 展 黄 定 了 
基础 。 面 向 传统 科学 计算 归纳 的 各 个 主题 ， 曾 为 应 用 
于 不 同学 科 的 数值 方法 (以 及 并 行 计算 ) 的 研发 工作 
提供 了 明确 的 路 线 图 ; 面向 模拟 智能 的 各 个 主题 同样 
不 局 限于 狭义 的 性 能 或 程序 代码 ， 而 是 激励 算法 、 编 
数据 结构 和 硬件 方面 的 创新 "1。 
3.2 控制 层 

控制 层 承 上 启 下 ， 在 整个 计算 系统 中 起 到 连接 和 
控制 算法 映射 与 硬件 执行 的 关键 作用 ， 在 现代 计算 机 
系统 中 表现 为 系统 软件 栈 。 本 文 仅 讨论 和 科学 模拟 相 
关 的 关键 组 件 。 模 拟 智 能 计算 系统 的 控制 层 的 变化 主 
要 来 自 2 个 方面 : CD 数值 计算 、 大 数据 和 人 工 智 能 的 
RIS, @ 底层 硬件 技术 可 能 发 生 的 颠 履 性 变化 。 近 
年 来 ， 由 于 科学 大 数据 的 急剧 增加 ， 在 科学 模拟 的 数 
值 计算 阶 段 ， 大 数据 软件 栈 逐 渐 被 超 算 系统 领域 所 关 
注 ， 只 是 相对 于 传统 的 数值 计算 ， 大 数据 软件 栈 是 完 
全 独立 的 ， 在 模拟 流程 上 属于 不 同 的 步 又。 因此 ， 基 
于 2 套 系统 的 软件 栈 是 基本 可 行 的 。 而 在 模拟 智能 阶 
段 ， 情 况 产 生 了 根本 上 的 变化 。 根 据 前 文中 所 表示 的 
问题 解法 描述 公式 产 FUooD,AC))， 人 工 智能 和 大 数据 
部 分 都 是 向 入 在 数值 计算 内 的 ， 这 种 结合 是 一 个 紧 籽 
合 的 模拟 过 程 ， 自 然 需 要 一 个 异 质 融 合 的 系统 软件 
栈 。 以 DeePMDm*I 为 例 ， 该 模型 包含 平移 不 变性 的 
悉 入 网 络 、 对 称 性 保持 操作 和 拟 合 网 络 3 个 模块 。 鉴 
于 体系 的 能 量 、 受 力 等 属性 不 以 人 为 定义 改变 〈 例 
如 ， 便 于 测量 或 描述 而 赋予 体系 中 各 个 原子 的 坐标 )， 
接 和 人 拟 合 网 络 进行 原子 能 量 和 受 力 的 拟 合 ， 就 能 得 到 


HH 


程 语言 、 


较 高 精度 的 拟 合 结果 。 再 考虑 模型 的 训练 数据 强 依赖 
于 第 一 性 原理 计算 ， 整 个 流程 是 一 个 数值 计算 和 深度 
学 习 紧 而 合 的 过 程 。 

因此 ， 系 统 软件 在 代码 生成 和 运行 时 执行 过 程 
中 ， 将 不 再 区 分 共性 核 函数 的 来 源 ， 即 不 再 区 分 是 否 
由 传统 人 工 智 能 、 传 统 数值 计算 或 根据 特定 问题 进行 
人 工 定制 扩展 得 来 。 相 应 的 ， 系 统 软件 一 方面 需要 针 
对 3 类 不 同 来 源 的 共性 核 函数 ， 提 供 易 于 扩展 和 开发 
的 编程 接口 。 另 一 方面 则 需要 对 这 3 类 函数 ， 在 代码 
编译 方面 和 运行 时 资源 管理 方面 ， 兼 顾 并 行 效率 和 访 
存 局 部 性 等 性 能 保障 ; 在 面向 不 同 粒度 的 计算 任务 
时 ， 能 够 逐 层 进行 融合 和 协同 优化 ， 发 挥 不 同类 型 体 
系 结构 处 理 器 的 最 佳 性 能 。 
3.3 处 理 层 

纵 观 数值 计算 阶段 到 模拟 智能 阶段 ， 一 个 驱动 技 
术 发 展 的 重要 因素 是 当前 硬件 技术 无 法 满足 计算 需 
求 。 因 此 ， 处 理 层 设计 首要 问题 是 : 表示 层 的 变化 
(如 对 称 性 、 基 元 ) 会 产生 全 新 的 硬件 体系 架构 吗 ? 
它们 是 基于 传统 专用 集成 电路 (ASIC) 实现 ， 还 是 超 
越 互 补 金属 氧化 物 半 导体 (CMOS) 一 一 从 高 性 能 计 
算 的 发 展 路 线 图 来 看 ， 这 也 是 未 来 Z 级 超 算 的 硬件 设 
计 要 考虑 的 核心 问题 。 可 以 大 胆 预测 ， 在 2035 年 左 
右 ，Z 级 超 算 可 能 会 出 现 。 尽 管 基于 性 能 和 可 靠 性 因 
素 的 考虑 ， 那 时 CMOS 平 台 仍 将 占据 主流 ， 但 一 些 核 
心 组 件 将 是 建立 在 非 CMOS 工艺 上 的 专用 硬件 。 

摩尔 定律 虽然 放 组 但 依然 有 效 ， 要 重点 解决 的 关 
键 难题 是 如 何 副 近 摩 尔 定律 的 极限 。 换 句 话 说 ， 如 何 
通过 软 人 硬件 协同 设计 的 手段 ,将 基于 CMOS 的 硬件 潜 
力 充分 挖掘 出 来 。 因 为 ， 即 使 在 性 能 优先 级 最 高 的 超 
算 领 域 ， 多数 算法 负载 所 获得 的 实际 性 能 仅仅 只 是 硬 
件 裸 性 能 的 极 小 部 分 。 回 顾 超 算 领 域 早期 发 展 阶段 ， 
其 基本 设计 哲学 就 是 软 硬 件 协同 。 未 来 十 几 年 ， 微 处 
理 需 迅速 发 展 的 “红利 ”将 耗 尽 ， 面 向 模拟 智能 的 计 
算 系 统 硬 件 架 构 应 该 回归 到 从 头 设计 的 软 硬 件 协同 技 
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术 上 。 一 个 突出 的 合子 是 如 前 所 述 的 分 子 动力 学 模 的 科学 挑 成 。 

拟 ，Anton 系列 5 是 一 个 从 零 开始 设计 的 超级 计算 机 

家 族 ， 可 以 满足 大 规模 长 时 间 尺 度 的 分 子 动力 学 模拟 

计算 ， 而 这 也 恰恰 是 对 微观 世界 探索 的 必要 条 件 之 1 Weaver W. Science and complexity. American Scientist, 
一 。 然 而 ， 最 新 的 Anton 计 算 也 只 能 对 基于 经 典 力 场 LL 0 Oe 

模型 实现 20 微 秒 的 模拟 ， 无 法 进行 第 一 性 原理 精度 的 2 Winsberg E. Computer simulations in science// Zalta E N, 
长 时 间 尺 度 模 拟 ; 然而 ， 后 者 才能 满足 多 数 实际 应 用 
(如 药物 设计 等 ) 需求 。 


参考 文献 


Nodelman U, eds. The Stanford Encyclopedia of Philosophy 
(Winter 2022 Edition). (2019-09-26) [2024-01-07]. https:// 


plato. stanford. edu/archives/win2022/entries/simulations- 


最 近 ， 作 为 模拟 智能 的 典型 应 用 ，DeePMD 模型 Science. 
在 传统 大 规模 并 行 系统 上 的 突破 证 明了 其 巨大 的 潜 3 Lynch S M. Using Statistics in Social Research: A Concise 


力 。 中 国 科 学 院 计算 技 术 研 究 所 超 算 团 队 ， 已 实现 了 Approach. New York: Springer, New York, 2013. 
170 个 原子 的 第 一 性 原理 精度 分 子 动力 学 的 纳 秒 级 模 4 Schelling T C. Dynamic models of segregation. The Journal 
m JA — TEETH Ic YZ 
HOU. ER. Fs e o SE OR E EU DS of Mathematical Sociology, 1971, 1(2): 143-186. 
ME. SE BR pj SAE FGA HS 人 
的 可 扩展 性 ,需要 在 运算 逻辑 和 通信 操作 上 有 极 钳 的 scientific discovery. npj Systems Biology and Applications, 
创新 。 本 文 认为 有 2 类 技术 可 以 预期 能 够 发 挥 关 键 作 2021, 7: 29. 
Hi: (D 存 算 一 体 架 构 通过 降低 数据 移动 的 延迟 来 提 6 Klahr D, Kotovsky K. Complex Information Processing: The 
运算 效率 ; (OO 硅 光 互 连 技术 ， 可 以 在 高 能 效 下 提供 Impact of Herbert A. Simon. New Jersey: Lawrence Erlbaum 

大 带宽 的 通信 能 力 ， 有 助 于 提高 并 行 性 和 数据 规模 。 Associates 989: 373-398. 
进而 ， 随 着 对 模拟 智能 应 用 广泛 而 深入 的 研究 ， 相 信 
未 来 将 逐步 形成 科学 模拟 领域 的 “新 浮 点 ”运算 单元 

A 

YY 


SE 02)[2024-01-07]. https://netlib.org/benchmark/hpl/. 


Nn 


Kitano H. Nobel Turing Challenge: Creating the engine for 


N 


Petitet A, Whaley R C, Dongarra J, et al. HPL—A Portable 


Implementation of the High-Performance Linpack 


Benchmark for Distributed-Memory Computers. (2018-12- 


8 Dongarra J, Heroux M A. Toward a New Metric for Ranking 
4 总 结 


High Performance Computing Systems. East Ave, 
ear sf aie ue S didus Livermore: Sandia National Laboratories, 2013. 
本 文 认为 ， 在 科学 模拟 的 现 阶段 ， 尚 处 于 模拟 智 
— n" N n 9 Dirac P A M. Quantum mechanics of many-electron systems. 
能 阶段 的 早期 ， 此 时 对 模拟 智能 的 使 能 技术 展开 研究 l | | 
Proceedings of the Royal Society of London Series A, 1929, 
至 关 重 要 。 在 一 般 科 学 研究 中 ， 独 立 的 概念 、 关 系 和 3 10733 


行为 可 能 是 易 理 解 的 。 但 是 ， 它 们 的 组 合 行为 会 导致 10 Lavin A, Krakauer D, Zenil H, et al. Simulation intelligence: 


不 可 预测 的 结果 。 深 入 了 解 复 杂 系 统 的 动态 行为 ， 对 Towards a new generation of scientific methods. (2021-11- 
于 许多 处 理 复杂 挑战 性 领域 的 研究 人 员 来 说 是 非常 宝 27)[2024-01-07]. https://arxiv.org/abs/2112.03235. 


贵 的 。 在 面向 模拟 智能 的 计算 系统 设计 中 ， 一 个 必 不 。 “1! RAE PEF MAh, 等 .科学 智能 (AI4S) ERREN 
可 少 的 环节 是 路 学 科 合作 ， 即 领域 科学 、 数 学 、 计算 。 URE ORE TAREA TA EAA, MIS 

Zhang L F, Li X Y, EWN, et al. Observation and Outlook on 
机 科学 与 工程 、 建 模 与 仿真 等 学 科 工 作者 之 间 的 协 the Global Development of AI for Science (AI4S). Beijing: 
作 。 这 种 跨 学 科 合作 会 构建 更 优 的 模拟 计算 系统 ， 形 AI for Science Institute, Beijing, 2023. (in Chinese) 
成 更 全 面 和 整体 的 方法 ， 去 解决 更 加 复杂 的 现实 世界 12 Zhang X, Wang L M, Helwig J, et al. Artificial intelligence 


24 |2024 年 .第 39 卷 .第 1 期 


面向 模拟 智能 的 计算 系统 | Sepe 


for science in quantum, atomistic, and continuum systems. 
(2023-07-17)[2024-01-07]. http://arxiv.org/abs/2307.08423. 

13 Voit E O. Perspective: Dimensions of the scientific method. 
PLoS Computational Biology, 2019, 15(9): e1007279. 

14Kuhn T S, Hawkins D. The structure of scientific 
revolutions. American Journal of Physics, 1962, 31(7): 
554-555. 

15 Wah B W, Li G J. A survey on the design of multiprocessing 
systems for artificial intelligence applications. IEEE 
Transactions on Systems, Man, and Cybernetics, 1989, 19 
(4): 667-692. 

16 Colella P. Defining Software Requirements for Scientific 
Computing, 2004, (2012-05-19) [2024-01-07]. https://www. 
slideserve.com/kiral/software-requirements. 

17 Asanovié K, Bodík R, Catanzaro B, et al. The landscape of 
parallel computing research: A view from berkeley. Technical 
Report No. UCB/EECS-2006-183. 


Berkeley: EECS 


Department, University of California, 2006. 


18 Zhang L F, Han J Q, Wang H, et al. Deep potential molecular 
dynamics: A scalable model with the accuracy of quantum 
mechanics. Physical Review Letters, 2018, 120(14): 143001. 

19 Jia W L, Wang H, Chen M H, et al. Pushing the limit of 
molecular dynamics with ab initio accuracy to 100 million 
atoms with machine learning// SC20: International 
Conference for High Performance Computing, Networking, 
Storage and Analysis. Atlanta: IEEE, 2020: 1-14. 

20 Shaw D E, Adams P J, Azaria A, et al. Anton 3: Twenty 
microseconds of molecular dynamics simulation before 
lunch// SC21: International Conference for High 
Performance Computing, Networking, Storage and Analysis. 
St. Louis: IEEE, 2021: 1-11. 

21 Guo Z Q, Lu D H, Yan Y J, et al. Extending the limit of 
molecular dynamics with ab initio accuracy to 10 billion 

27th ACM SIGPLAN 


atoms// Proceedings of the 


Symposium on Principles and Practice of Parallel 


Programming. Seoul: ACM, 2022: 205-218. 


© 4 £2) 院 刊 | 25 


| 专题 : 大 力 推进 科研 范式 变革 


Computing system for simulation intelligence 


TAN Guangming’ JIA Weile WANG Zhan YUAN Guojun SHAO En SUN Ninghui’ 

(Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, China ) 
Abstract This study refers computer simulation in scientific research to as scientific simulation. Based on its narrow and broad 
definitions, this study divides scientific simulation into three stages: numerical computation, simulation intelligence, and science brain. 
Now, scientific simulation is entering the era of simulation intelligence, i.e., driven by scientific big data and artificial intelligence 
technology, scientific simulation is shifting from traditional numerical simulation to simulation integrated with artificial intelligence. In 
order to understand what the right computing system for simulation intelligence is, the design guidelines, basic methods, and key 
technical problems are discussed. 


Keywords scientific simulation, simulation intelligence, artificial intelligence, computing system, Zetta-scale computing 
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